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AMPLITUDE AND PHASE FLUCTUATIONS FOR GRAVITATIONAL WAVES PROPAGATING THROUGH 
INHOMOGENEOUS MASS DISTRIBUTION IN THE UNIVERSE 
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ABSTRACT 

When a gravitational wave (GW) from a distant source propagates through the universe, its am- 
plitude and phase change due to gravitational lensing by the inhomogeneous mass distribution. We 
derive the amplitude and phase fluctuations, and calculate these variances in the limit of a weak 
gravitational field of density perturbation. If the scale of the perturbation is smaller than the Fresnel 
scale ~ 100pc(//mHz) -1 / 2 (/ is the GW frequency), the GW is not magnified due to the diffraction 
effect. The rms amplitude fluctuation is 1 — 10% for / > 10~ 10 Hz, but it is reduced less than 5% 
for a very low frequency of / < 1CU 12 Hz. The rms phase fluctuation in the chirp signal is ~ 1CU 3 
radian at LISA frequency band (10 -5 — 1CU 1 Hz). Measurements of these fluctuations will provide 
information about the matter power spectrum on the Fresnel scale ~ 100 pc. 
Subject headings: gravitational lensing - gravitational waves - large-scale structure of universe 



1. INTRODUCTION 

The inspiral and merger of SMBHs (Super Massive 
Black Holes : 10 4 — 10 7 M Q ) is one of the most promis- 
ing candidates for LISA (Laser Interferometer Space An- 
tenna), which will be launched around 2014. This detec- 
tor is sensitive in a frequency range of 10 —10 Hz and 
can measure SMBH mergers at cosmological distances 
with a high signal-to-noise ratio. The binary SMBH sys- 
tems have recently been called "cosmological standard 
sirens" (Holz & Hughes 2005). This is because the dis- 
tance to the source r s can be directly measured by using 
the relation A oc / /(/ 3 r s ), where / is the GW frequency, 
/ is its time derivative and A is the amplitude at the 
detector (Schutz 1986). The distance can be determined 
with less than 1% accuracy if the direction is determined 
by identifying its electromagnetic counterpart (e.g. Cut- 
ler 1998; Seto 2002; Hughes 2002; Vecchio 2004). 

But in practice, the distance cannot be determined 
with such high accuracy because of the gravitational lens- 
ing caused by inhomogeneous mass distribution in the 
Universe. Recently, Holz & Hughes (2005) and Kocsis 
et. al. (2005) discussed the effects of lensing magnifi- 
cation (or demagnification) on determining the distance 
to SMBH mergers. They concluded that lensing errors 
are 5 — 10%, which is greater than the intrinsic distance 
error. 

In weak gravitational lensing, the magnification is 
fi ~ 1 + 2k, where k is the convergence. The rms conver- 
gence fluctuation was derived by Blandford et al. (1991), 
Miralda-Escude (1991) and Kaiser (1992) on the basis of 
linear perturbation theory. In this previous work, the 
geometrical optics was assumed. 

But for the lensing of GWs, since the wavelength is 
much longer than that of light, the geometrical optics is 
not valid in some cases. Recently, suggested by Macquart 
(2004) and Takahashi, Suyama & Michikoshi (2005), if 
the scale of the density perturbation k~ 1 is smaller than 
the Fresnel scale rp ~ (Xts) 1 ! 2 (A is the wavelength), the 
wave effects become important. This condition is rewrit- 
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ten as k' 1 < 100pc(//mHz)- 1 / 2 (r s /10Gpc) 1 / 2 . In such 
a case, the incident wave does not experience perturba- 
tion and its amplitude is not magnified. 

In this paper, we consider a situation in which GWs 
propagate through the density perturbation of CDM 
(cold dark matter) and baryon. In the wave optics, the 
lensing affects not only the amplitude but also the phase, 
hence we discuss the lensing effects on both of them. We 
use units of c = G = 1. 

2. GRAVITATION ALLY LENSED WAVEFORM 

The background metric is the Friedmann-Robertson- 
Walker (FRW) model with a gravitational potential of 
f/(<C 1). The perturbed FRW metric for a flat uni- 
verse is given by 

dS 2 = g^ v dx"dx v = a 2 {rj)gB v dxi x dx v 



■■a 2 (v) [- (1 + 2U) dr] 2 + {1 - 2U) dx 2 



(1) 



where r\ is the conformal time, the scale factor is normal- 
ized such that a = 1 at present, and gr[ v is the perturbed 
Minkowski metric. The line element is dx 2 = dr 2 +dxj_, 
where r is a radial coordinate (a comoving distance) from 
the observer, r{z) — J Q dz'/H(z'), while xj_ is a two- 
dimensional vector perpendicular to the line-of-sight. We 
show the lens geometry for the observer and the source 
in FigQ The observer is in the origin of the coordi- 
nate axes, while the source position is x s = (r a , Xj_) with 
|x^|<r s . 

Since the propagation equation of GW is (1) confor- 
mally invariant if the wavelength is much smaller than 
the Hubble radius (see Appendix) and (2) the same as the 
scalar field wave equation (Peters 1974), we use the scalar 
field <t> propagating under the Minkowski background 

= or 



9, 



p V . The basic equation is <9 M (■s/—g B g l gd v 

(v 2 +lu 2 ) $ = Auj 2 U4>, 

where u is the frequency at the observer 1 and 
the Fourier transform of <j>(r),x). 



(2) 



(id, x) is 



1 w is the same as the comoving frequency u c , since u) c = au 
and a = 1 at present. 
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Fig. 1. — Gravitational lens geometry for the observer and the 
source, r is the comoving distance from the observer, r(z) = 
Jq dz' I H(z'), and is a two-dimensional vector perpendicular 
to the line-of-sight. The source position is x a = (r B ,x s ± ) with 
<S r s . (f>° is the incident wave, and <jp + </> 1 is the lensed wave. 



We take 0° as the the incident wave emitted from the 
source, which is the solution of Eq.((2J) in the unlensed 
case U = 0. We use the spherical wave as 0° : <^> (u>, x) oc 
e iw|x-x s |^| x _ Xs | Including the effect of U on the first 
order (Born approximation), the gravitationally lensed 
wave at the observer is given by ijTakahashi et al. 2005) 



with 



4 obsM = 

7T 



e M x l 

z^C/(x)/(w,x). 
x 



(3) 



(4) 



where (p 1 represents the effect of lensing <C |^°|)- 

The incident wave 4>° is gravitationally lensed at x — 
(r,xx) and changed into the lensed wave + cf) 1 (see 
FigGJ. 

Let us define K and S as (Ishimaru 1978, Ch.17), 



K(u) = Re 



(5) 



S(u>) = Im 

"ob B MJ L%bs( W ) 

Then we have 

^M = [l + K( L u)]4> o Me iS ^\ (6) 

from Eq.(J3J). Hence K means the magnification (K > 0) 
or demagnification (K < 0) of the wave amplitude, while 
S means the phase shift due to the lensing. Hereafter, 
we call K amplitude fluctuation and S phase fluctuation, 
respectively. 

Using the Fourier transform of the potential, U(r, k) = 
/ d 3 xU(x.)e lkx , with |xj_| <C r and |x^| <C r s , the result 
in Eq.(0J is reduced to 



^bs( 



-2iw 



dr 



(27T) 



U(r,k) 



x exp 



— ik r r — i — kj 



.r(r s 



2u)r. 



(7) 



where k r and kj^ are the radial and perpendicular compo- 
nents of k. In particular for high frequency limit u> — > oo, 
K and S are rewritten as 



dr ■ 



-ViC/(r, -xi), 



S{uj)^Lut d = -2lu drU(r, — xl). 



(8) 



Here, k is the convergence field along the line-of-sight to 
the source and td is the gravitational time-delay. The 
above results are consistent with that in wea k gravita- 
tional lensing (Bartelmann & Schneider 2001). 

3. AMPLITUDE FLUCTUATION 

3.1. Variance in the Amplitude Fluctuation 

In this section, we derive the variance in the ampli- 
tude fluctuation K. The gravitational potential satisfies 
Poisson's equation (Peebles 1980) 



U(r,k) 



2a(r) 



(9) 



where 6 is the density perturbation. The fluctua- 
tion of S is characterized by the power spectrum : 
(5(r,k)6*(r',k')) = (27r) 4 P*(r, k)S 3 (k - k')S{r - r'). 
Then, the correlation in the potentials U (r, k) and 
U(r',k') is 



(U(r,k)U(r',k')) = 



2a(r) 



x (2tt) 4 P s (r, fc)<5 3 (k - k')6(r - r'). 



(10) 



To calculate Ps, we use the linear power spectrum 
ijEisenstein fc Hu 19 99) with the nonlinear correction 
of Peacock & Dodds (1996). We adopt a COBE- 
normalized 2 , scale invariant (n — 1) power spectrum 
in a flat ACDM cosmology with fl b = 0.04, fl Q = 0.3, 
A = 0.7, and H = 70 km s _1 Mpc _1 . 

The variance in the amplitude fluctuation is given from 
Eqs. ©, and (JUJ as 



An 



o a in 



2 / V 

r z (r s - r) 



d 2 k±P S {r,k±)F K (u J ,r,k 1 _) (11) 



where the filter function Fk is defined as 



K 



sin (4fcj/2) 
r 2 F k\/2 



r (r s 



(12) 



Here, rp is called the Fresnel scale ijMacauart 2004'). 
This is roughly given by 



120pc 



( / 



V mHz 



-1/2 



r (r s -r) /r s 



lOGpc 



1/2 



(13) 



where / = lo/2tt. We show the filter function Fk as a 
function of rpk± in FigJH From this figure, Fk = 1 if 
the scale of the density perturbation ~ kj 1 is larger than 

2 It corresponds to erg = 0.8 which is the present amplitude of 
the mass fluctuation at 8h~ 1 Mpc. 
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Fig. 2. — Filter function Fk as a function of rpk. 



the Fresnel scale rp, while Fx decreases rapidly in pro- 
portion to (fi?fcj_)~ 4 if kj 1 -C rp. Hence, the amplitude 
fluctuation is affected by a density perturbation of scale 
larger than the Fresnel scale. In the geometrical optics 
limit Fx — + 1, the result is the same as (k 2 ). 

The reason why a critical scale is rp can be ex- 
plained as follows : For the lensing by a compact ob- 
ject with a mass M, if the wavelength is larger than the 
Schwarzschild radius, A > M, the diffraction effect be- 
comes important (see Takahashi & Nakamura 2003, and 
references therein). Inserting this condition to the Ein- 
stein radius rp = [4M(r s — r)/?^] 1 / 2 , we have the Fresnel 
scale in Ea. (fT2*|l . 

3.2. Results 

We show the rms (root-mean-square) amplitude fluc- 
tuation Ax in Ea. (|f I \ as a function of the frequency in 
Fig El The source redshift is z s = 3. The solid line is 
Ax and the dashed line is the geometrical optics limit 

1/2 

^k 2 ) . The difference between these two lines is small 
for the LISA frequency band (10~ 5 to 10 _1 Hz). But for 
very low frequency / < 10 _ Hz, Ax is clearly smaller 
than (k 2 ) 1 / 2 . Thus the frequency band for the pulsar 
timing array / « 10 -9 Hz (e.g. Jenet et al. 2005), Ax 
should be used instead of (k 2 ) 1 / 2 . The wave amplitude is 
magnified due to lensing by the density perturbation on 
a scale larger than the Fresnel scale (see Sec.3.1). The in- 
tegration of Eq. (|f I \ is mainly contributed on the scale of 
0.1 - 1 Mpc (since Psk 2 has a peak there). For / < 10~ 10 
Hz, the Fresnel scale is r F > 0.4Mpc(//10 _10 Hz) _1 / 2 
which is larger than the contributed scale (0.1 — 1 Mpc), 
and hence the amplitude is not magnified. 

FigQJis the same as FigO but as a function of z s . The 
dashed line is the geometrical optics limit (k 2 ) 1 / 2 , the 
solid (dotted) line is for the frequency of 10~ 10 (10 -12 ) 
Hz. For z s = 1 — 10 the rms amplitude fluctuation is 
1 — 10% for / > 10~ 10 Hz, while it decreases less than 
5% for / < 10~ 12 Hz. 




frequency (Hz) 



Fig. 3. — The rms amplitude fluctuation Ajf as a function of the 
frequency. The source redshift is z B = 3. The solid line is A^, and 
the dashed line is the result in the geometrical optics limit (ft 2 ) 1 / 2 . 
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Fig. 4. — Same as Fig[3] but as a function of z B . The dashed line 
is {k 2 ) 1 / 2 , and the solid (dotted) line is Ape for the frequency of 
10~ 10 (1CT 12 ) Hz. 



3.3. Diffraction effect in the amplitude fluctuation 

To investigate the diffraction effect in the amplitude 
fluctuation, we expand K(lu) in terms of 1/w. From 
Eqs.© and Q, we have 

K(u) = 2lj j " dr sin Qr^V?^ U(r, ^x s ± ), (14) 

where sin(r 2 ,V 2 L /2) = ££° =1 (-l) n+1 (r|V 2 L /2) 2 "- 1 
/(2n — 1)!. The leading term (n = 1) is the convergence 
K in Eq.©. The n-th term (n > 2), being of the order of 
(rfrV^) 2 ^™" 1 ^ ~ (rp/k^) 41 - 71 ^ 1 ^ k, is a correction term 
due to the diffraction effect. The perturbation on a scale 
smaller than the Fresnel scale affects the diffraction. 

4. PHASE FLUCTUATION 
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4.1. Effects of the Phase Fluctuation in Chirp Signal 

Next, we discuss the effects of the phase fluctuation. 
We consider the inspiraling BH binaries as the sources. 
As the binary system loses its energy due to gravitational 
radiation, the orbital separation decreases and the orbital 
frequency increases. Hence, the frequency of the gravi- 
tational waves increases with time (df/dt > 0). This 
is called a chirp signal. We consider the frequency to 
be swept from f\ to f 2 . For the binary masses M\ and 
M 2 at redshift z s , the frequency of lyr before the final 
merging is 



h = 4.1 x 10 



-5/ M t 



-5/8 



Hz 



(15) 



V io 6 m q 

where M z = {M 1 M 2 f /5 {M l + M 2 )- l l h (l + z s ) is the 
rcdshiftcd chirp mass. The frequency at the final merging 
is 

/ M z 



h = 4.4 x 10" 



V 10 6 M Q 



Hz 



(16) 



where M z = (Mi + M 2 )(l + z s ) is the redshifted total 
mass. The difference in the phase fluctuation between f\ 
and f2 in the chirp signal is important. 

The phase fluctuation S is reduced to the time delay 
utd in the geometrical optics limit from Eq.©. But, 
we note that the time delay is physically unimportant, 
since it means an arrival time shift and hence it does not 
change the waveform. Hence we use S — totd instead of 
S as the phase fluctuation. This quantity S — totd has 
larger value for smaller frequency. We define As as the 
rms phase difference between the two frequencies, u>\ and 
u 2 ■ 

A|(wi, u> 2 ) = ([(S(ui) - ujit d ) - (S(lu 2 ) - u 2 t d )] 2 ). 

(17) 

This is the same as A 2 ^ in Eq. 1)11(1 but the filter function 
is replaced with 



(4i*l/2) 



F2 



kl/2) 



?1 k ± /2 



F2 



kl/2 



(18) 

where rpi and rp 2 are the Fresnel scales in Eq. <|12|) for 
w\ and u> 2 , respectively. 

We show the behavior of first and second terms in 
Eq.fp, I [cos(rf,fc 2 /2) - l] /(r 2 F k 2 /2)\, as a function of 
rpk in Fig|3] From this figure, the function peaks at 
rpk ~ 1. Hence the phase fluctuation S — Lutd is affected 
by the density perturbation of the Fresnel scale. In the 
chirp signal, as the frequency increases, the GW feels the 
perturbation of the smaller scale. 

4.2. Results 

In table ^ we show the rms phase differences As in 
Eq.jnjl for the LISA frequency band (10~ 5 to 10" 1 Hz) 
with z s — 1,3 and 10. We consider the frequency to be 
swept from f\ to f 2 in the chirp signal. The values arc in 
units of radian. This table shows that the typical values 
of As are w 10~ 3 radian. The results weakly depend on 
fi if fx *C fi- This is because S — utd in Eq. l(T7|) has 
larger (smaller) value for smaller (higher) frequency. 

In FigEJ As is shown as a function of f\ in the limit 
of f 2 — > 00. The source redshifts are z s = 1 (dotted 
line), 3 (solid) and 10 (dashed). In order to study the 




Fig. 5.— I [cos(r|,fc 2 /2) - l] /(r 2 F k 2 /2)\ as a function of r F k. 



TABLE 1 

RMS PHASE DIFFERENCES IN CHIRP SIGNAL f 1 — * fa Hz WITH 

2 S = 1, 3, 10. The unit is radian. 



h 




h (Hz) 
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behavior of As, we assume a single power law for the 
power spectrum, P(r,k) oc k n . The index is n ~ —2.7 
at fc _1 ~ 100 pc. Inserting this P(r,k) into Eas. (fTT|l 



and (HHJ, we have As oc cj 1 ™ + ]l oc cj^ ' 18 . With this 
argument and the results in table ^ As in f\ -C f 2 is 
roughly fitted by 

-0.18 



3 x 10- 



rad 



fi 



10- 4 Hz 



for z s = 3. (19) 



The above value, 3 x 10 3 rad, is replaced by 1(5) x 10 3 
rad for z s = 1(10). 

4.3. Implications for GW Observations 

In the matched filtering analysis, the phase of wave- 
form can be measured to an accuracy approximately 
equal to the inverse of the signal-to-noise ratio, (S/N)" 1 . 
Here, S/N is typically « 10 3 for the SMBH mergers de- 
tected by LISA. But including the effect of the phase 
fluctuation, the phase cannot be determined with an ac- 
curacy of less than As- Hence if S/N is larger than 
A^ 1 = 10 3 (As/10 _3 rad), the phase fluctuation becomes 
important and the phase of the waveform can be deter- 
mined with an accuracy of ~ Ag 1 . 

4.4. Diffraction effect in the phase fluctuation 
Similar to Sec. 3.3, we rewrite the phase fluctuation as 

1 
2 



S(w) = -2lu 



dr cos I ^r|V 2 



± , U(r, -xl). (20) 
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Fig. 6. — As as a function of /i in the limit of /2 — » oo. The 
source redshifts are z s = 1 (dotted line), 3 (solid) and 10 (dashed). 



through the strong density fluctuation or pass near mas- 
sive compact objects (i.e. strong lensing) . Hence we have 
some comments about the validity of the Born approxi- 
mation. 

For the amplitude fluctuation K, the maximum of 
K is the convergence k from FigEI For the thin lens 
plane at ri, k is the surface density of the lens S di- 
vided by the critical density Scr = (l/47r)r s /r^r^ s ~ 
2 x l0 3 M Q pc- 2 [{r s /r e re s )/Gpc~ 1 ] (here r £s = r s - r e ). 
Hence if the GWs pass through high density region 
£ > Ecr such as a core of galaxy or cluster, the Born ap- 
proximation breaks down and one should use the Kirch- 
hoff diffraction integral to obtain the exact waveform 
(Schneider, Ehlers & Falco 1992, Sec.4.7 and 7). 

If the gravitational potential is a Gaussian random 
field, one can obtain exact solutions of the correlation 
functions of the lensed waveform (Macquart 2004; see 
also Ishimaru 1978, Ch.20). If there are many compact 
lens objects and the GWs are scattered several times, 
the exact lensed waveform can be obtained by using the 
multiple lens-plane theory in the wave optics (Yamamoto 
2003). 



In u — > 00, the above result is wtd in Eq.jnj). To study 
the behavior of the correction terms, we consider a simple 
model. We assume the density fluctuation is confined in 
a thin lens plane at distance n and is fitted by a single 
power law near the Fresnel scale : U (r, k) cx k°[8(k r )8(r— 
ri) where a is an index. Then, from Eqs.(|5jl and (JJJ, we 
have 

S(u) =ujt d + Cuj a/2+2 . (21) 

The coefficient C is constant and the second term is of the 
order of 10 -3 rad at the LISA frequency band. Hence, 
by measuring the phase fluctuation directly, one would 
obtain the constant C and the index a. 

5. VALIDITY OF THE BORN APPROXIMATION 

Throughout this paper we assume the weak density 
fluctuation and employ the Born approximation to dis- 
cuss the lensing effects on the waveform. Since the vari- 
ances of the amplitude and the phase fluctuations are 
much smaller than one, this approximation is valid in 
almost all cases. But in a few cases, the GW may pass 



6. CONCLUSION & DISCUSSION 

We have discussed the lensing effects on the amplitude 
and phase of a waveform. The rms amplitude fluctua- 
tion is 1 — 10% for / > 10~ 10 Hz, which is the same as 
the result in weak lensing. But for a very low frequency 
of / < 10~ 12 Hz, it decreases to less than 5%. In the 
chirp signal, the phase fluctuation is typically 10~ 3 ra- 
dian at the LISA frequency band. The phase cannot be 
measured with an accuracy less than this value. 

The power spectrum P$ has been measured down to a 
small scale ~ several x 0.1 Mpc from the Lya forest 
(e.g. Zaroubi et al. 2005). In this paper, we assume the 
formula of P$ is valid down to the Fresnel scale fc _1 ~ 100 
pc. If the amplitude or the phase fluctuation is measured 
in future, the constraints for P$ at ~ 100 pc could be 
obtained. 



I would like to thank Naoshi Sugiyama and Naoki Seto 
for useful comments and discussion. I also thank anony- 
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electromagnetic scattering. 



APPENDIX 
CONFORMAL TRANSFORMATION 
The relation between the Einstein tensor G^ v for the metric g M „ and G^ v for g^ v {= a 2 g^ v ) is given by l)Wald 1984|) . 

G MV = G^y - IVySv In a + 2 (V M In a) (V„ In a) 

+g„ v (V CT In a) (V CT In a) + 2 ffjtl!/ V <T V (T In a. (Al) 

Let us consider the linear perturbation h llv (= a^h^) in the background metric g^ v {= o?g^ v ) given in Eq.(|T]l: 

Say = 9% + V = « 2 (9%, + V) ■ ( A2 ) 
Inserting Ea. (|A2|l into (|A1|I . we obtain the linearizing Einstein equation 



<5G>„ = 5G^ U + 2—ST° 



a" a' 2 
2 2 



S (g^g™) - 2^-5 (W"T£ P ) , 



(A3) 



where a' = da/dr], T is the Christoffel symbol, and 5 means the perturbed component. The first term 5G^ V is of the 
order of \h^ v \/\ 2 ( A is the wavelength of the gravitational wave ). The other terms on the right-hand side are roughly 
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|/1/ik|/(AAh) or \h^ u \/)? H where \r is the horizon scale, since a' = a 2 H ~ a 2 Hence, if A <C \h the propagation 
equation for the gravitational wave is conformally invariant, i.e. 8G^ V = 8G^ V = 0. 
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